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Abstract 

The potentiality of a Giant Shower Array to low energy gamma rays from 
gamma ray bursts is discussed. Effective areas are calculated for different 
scenarios and the results are encouraging. If gamma ray bursts have a spec- 
trum which continues in the high energy gamma ray region, the Pierre Auger 
Observatory will be able to detect it. 



I. INTRODUCTION 

Gamma Ray Bursts (GRB) are one of the most intriguing mysteries in the Universe and 
models abound trying to explain them. A recent review of models is given in ref . [|l| . In ref . 
a compilation of more than 100 models is given. A complete phenomenological review 
with many references is given in ref |Q. 

Recently, thanks to the accurate measurement of the position of a GRB by the BeppoSAX 
satellite |Q it was possible to measure a source counterpart in other frequencies (X rays, 
optical and radio bands) and a measurement of the redshift [Q was given for the GRB970228 
(0.86 < z < 2.3) implying a truly cosmological distance. 

In this respect it will be most interesting to measure the flux of high energy photons from 
gamma rays burst. We should remember that one does not expect high energy photons if 
GRB are cosmological. Even a 100 GeV photon flux would suffer a considerable attenuation 
for cosmological distances (~ Gpc). Therefore, looking for high energy photons would be 
an invaluable tool in discriminating cosmological scenarios from extended halo or mixed 
models. 

In this letter we discuss the potential of a giant surface array, specifically the Pierre 
Auger Project |0 to measure bursts fluxes of low energy ( GeV - TeV) photons. Studies 
of this type have been done for a variety of detectors and limits have been put by different 
collaborations, see for instance [0,H. All present results are negative and only upper limits 
could be put on fluxes on high energy. An interesting result is, however, the 10 a candidate 
event detected by EAS-TOP 0. Although the event could not be associated to any GRB's 
this could be an indication of the " delayed phenomena" . In this work, the acceleration 
of the highest energy cosmic rays is related to GRB. Delayed high energy photons appear 
naturally as a consequence of the propagation and cascading of the cosmic rays through 
the photon (CMBR and IR) fog. The delay time depends on the intergalactic magnetic 
structure and could be of the order of hours to days (or even years). 

Recently, a work similar to ours was written by DuVernois and Beatty [|T0| were calcu- 



lation of the effective area for the Auger detector was done. However, our results indicate 
that their effective area is overestimated by a factor of 15 with respect to us. We have been 
unable to trace back the cause of discrepancy. 

II. SIMULATION 

We have run several thousands photon initiated showers for different energies using the 



program Aires |Tl|. Threshold energies for both muons and e.m. particles have been set 
above the respective threshold for Cerenkov light production in water. The ground depth 
was set to 850 gr/cm^ 0. In table | we can see some parameters of the run. The thinning 
level is always chosen so as to guarantee that all particles in the shower are kept until they 
reach the threshold value. 

For every shower we have done the following simulation. At ground level an infinite grid 
is built. Each cell in the grid has an area of a = 10 m^ equal to the area of the detectors 
in the Auger Project. For each shower we compute the number of particles that reach each 
cell at the ground. If the number of particles is above a given number, k, we will count it 
as a trigger. Results of the simulation are given in the table |. Let Nk{E) be the number of 
triggers for showers of energy E. Then a convenient parametrization is given by: 

NkiE) = No{k) E'-'', (1) 

where A^o(fc) depends on k but it is independent of the energy. In fig.Q we can see the result 
of our calculation of Nk for different values of fc as a function of energy. Also plot is the 
above parametrization. We can see that for medium energies { 10 < E < 100 GeV) the 
parametrization is good. At high energies it overestimates the values of A^^. 

Given the value of Nk we can calculate the effective area for detecting low energy showers 
as follows. A shower of some given energy will have an "effective area" given by: 

As{E) = aNk{E), (2) 



where a is the area of each detector. This expression for the effective area reflects the low 
energy character of the photon initiated shower which has shower maximum position much 
higher than the ground level, giving a surface distribution of particles bearing no relation 
to the initial shower axis direction. 

For an array of detectors, we have that the probability of a shower to trigger a detector 
is given by the ratio of the shower area by the inter detector area i.e.: 

P(E) = ^, (3) 

where / is the separation distance between detectors. So the total effective area will be: 

A^QiE) = At PiE), (4) 

where A^ ~ NdI"^ ~ 3000 km^ is the total area covered by the detector and N^ is the 
number of detectors. Finally we get: 

A^QiE) = aNnNk{E). (5) 

From fig. |1] we can estimate the effective area for the Auger Project. At i? ~ 100 GeV 
we have A q ~ a N^ ~ 1.6 10^ m^. Even at so low energies as -E ~ 10 GeV we have an 



effective area of A^^ ~ 16 m , i.e. bigger than the effective area of EGRET |12|. This result 



is in contradiction with the result in ref . [^ . In fig. [^ we show their result as a continuous 



line. We can see that their result is ~ 15 times higher than ours even in the most optimistic 
scenario of A; = 1. In comparison we also show the EAS-TOP result 0] which agrees with 
ours. This should be expected since the altitude for EAS-TOP is similar to the projected 
Pierre Auger (850 gr/cm^) altitude and the area of each detector is equal (10 m^). We 
conclude that our results are correct. 

Due to the large distance between detectors for the Auger Project (/ ~ 1.5 km) the array 
always have to operate in single counting mode. Even at the highest energies (~ 100 TeV) 
the shower will certainly trigger on at most one detector and no correlation on neighbour 
detectors should be expected. This is in contrast to smaller arrays where the small distance 



between detector allows to measure correlation between neighbour detectors see for instance 

i- 

For the angular dependence we have simulated 1000 showers at fixed energies and at 

different angles. In table ^ we show the result of the simulation. We can parametrize the 

results in the form: 

Nk{e) = No cos(^)", (6) 

with a ~ 9 which agrees with previous results ||TO|J^. We assume that this dependence 
is valid for other energies. With the result of our calculation we are able to calculate the 
effective area for GRB with arbitrary spectrum and for arbitrary arrival direction. 



TABLES 



Energy (GeV) 


Nsh 


Nk=i 


Nk=5 


1 


10^ 


3 


1 


5 


10^ 


40 


8 


10 


10^ 


250 


62 


20 


10^ 


767 


216 


25 


10^ 


1.12 10^ 


297 


50 


10^ 


4.45 10^ 


1.13 10^ 


100 


10^ 


1.52 10^ 


3.72 10^ 


150 


10^ 


2.89 10^ 


6.88 10^ 


200 


10^ 


4.79 10^ 


1.16 10^ 


250 


10^ 


6.91 10^ 


1.64 10^ 


300 


10^ 


9.36 10^ 


2.23 10^ 


500 


10^ 


2.02 10^ 


4.80 10^ 


Energy (TeV) 


Nsh 


Nk=i 


iVfc=5 


1 


10^ 


5.75 10^ 


1.39 10^ 


10 


5 103 


6.27 10^ 


1.98 10^ 


20 


10^ 


2.61 10^ 


9.53 10^ 


50 


102 


6.09 10^ 


2.67 10^ 


100 


102 


1.09 10^ 


5.35 10^ 


500 


70 


2.29 10^ 


1.37 10^ 



TABLE L Shower simulation parameters. 



Energy 


= 100 GeV 










e 


cos(^) 


Nsh 


Nk=i 


Nk=5 







1 


10^ 


15234 


3718 




10 


0.985 


103 


1295 


362 




20 


0.949 


10^ 


786 


196 




30 


0.866 


10^ 


334 


90 




40 


0.766 


103 


75 


19 



TABLE II. Angular dependence on trigger number. 



Let's assume a burst with a spectrum dN/dE = ^o{E/Eo) ^ between -E'max and -E'j-QJn 
occurs during a time T. Then the number of triggers in excess observed will be 

/■^max E ~'^ , , 

S = T dE%— A^Q{E,e). (7) 

J E ■ -C/n 

mm " 

In the same time the number of background triggers will be A^ = z/T, where u is the 
background trigger rate. And the statistical significance is given by: 

_ ^ ^ a<l>oTNDaj,9) 
'''' ^ VTVNE ^ ^ 

where C(7, 0) is given by 

Ch,e)= dEE~m,iE,e). (9) 

JE ■ 

mm 
Therefore a limit with a ricr confidence level can be obtained from no observation for fluxes 
bigger than 



For the Auger detector the single counting ratio will be about 2.5 kHz and the integration 
time can be of order 1 s, although given the uncertainty in the time profile for high energy 
photons, specially if the delayed phenomena is general, this time should be refined. Thus 
we get, assuming a spectrum index of 1.5 in the region of 1 GeV to 1 TeV: 

$0 = 5.1 X 10-V 2^^;, ,, - (11) 

cm^s GeV 

We should use Ua > 10 in order to avoid accidental triggers. With this values we would get 
upper limit fluxes which are competitive with actual Cerenkov detectors, see for instance 
ref. H. 

The present technique will not be able to give any indication of the arrival direction of 
these photons, contrary to Cerenkov detectors, therefore the actual detection will be possible 
only on the basis of timing considerations and correlations with other experiments. An 
advantage, however, is that it is sensitive to a much larger range of energies than Cerenkov 

8 



detectors: from lower energy photons to photons of hundreds of TeV and more (hmited 
only by the source), allowing, in principle, to fill the gap between EGRET and Cerenkov 
detectors. 

Acknoledgements 
The authors thank the IFUSP and IFT for providing us with computer resources. RAV 
thanks the IFT for its kind hospitality. This work was supported by CNPq (COE and PLS) 
and FAPESP (RAV). 



REFERENCES 



[1] M.J. Rees, |astro-ph/970TT62 . 



[2] R.J. Nemiroff, |astro-ph/940^im . 



[3] G.J. Fishman and C.A. Meegan, A.R.A. & A 33 (1995) 415. 

[4] E. Costa et al. (1997) lAU Circ. No. 6572. 

[5] M.R. Metzger, et al. (1997) lAU Circ. No. 6655. 

[6] See the Design Report for the Auger Project. Fermilab, Nov. 1996. See also Nucl. Phys. 
(Proc. Suppl.) 28B (1992), ed. M.Boratav, J. Cronin and A. A. Watson. 

[7] M. Aghetta et al. (EAS-TOP coll.) Astro.J. 469 (1996) 305. 

D.E. Alexandreas et al. (Cygnus coll.) Astro.J. 426 (1994) L1-L3. 

[9] E. Waxman and P. Coppi, IASSNS-AST-96-22, |astro-ph/9603T14 



[10] M.A. DuVernois and J.J. Beatty, Gamma-ray burst detection with the Auger 



Surface Array, DRAFT GAP-97-068, available in the WWW at |http://www^ 



auger. fnal.gov:82/GAP_Index.htnii 



[11] S.J. Sciutto, AIRES: A minimum document, DRAFT GAP-97-029. 
[12] G. Kanbach et al, Space Sci. Rev. 49 (1988) 69. 
[13] M.G. Baring, |astro-ph/97lT25^ . 



10 



FIGURES 



11 



-I — I — I I I I III 1 — I — I I I I III 1 — I — I I 1 1 III 1 — I — I I 1 1 III 1 — I — I I I I III 1 — I — I I I I III 1 — I — I I 1 1 III 



10 



3 _ 



B 
B ^ 

B ♦ 
B ^ 



iz; 



10 



_ 



10 



-3 _ 



-I I 



10 



-1 



m 




Dk=l 
O k=4 

+ k=5 



_l I I I I I I I I I I I I I I I L. 



10^ 10^ 



10^ 10^ 10^ 

E (GeV) 



10^ 10 



12 



FIG. 1. Number of triggers as a function of photon energy and k (see text) for photons arriving 
vertically. Also shown is the parametrization of DuVernois and Beatty (continuous line) and the 
parametrization of EASTOP (dashed line). 



13 



oa 
O 
v 

M 
2 



1.5 


■1 




...... 


,...., 


............. 


1.0 


■ 




D k=l 




/ : 




■ 




O k=4 




/ 




■ 




+ k=5 




/ 


0.5 


- 






^ 






■ 






, . - 


. - - ' i 


rv rv 


1 


D 

. ♦. ■ : 


■ : ' 1 . . . . 


----r" 

1 .... 1 


.... 1 .... 1 . . 



0.75 



0.8 



0.85 0.9 

cos(0) 



0.95 



14 



FIG. 2. Number of triggers as a function of arrival angle and k for photons of 100 GeV. Also 
shown are the cos{9)" fit. 
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